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Abstract 

We present a thermal inflation model that incorporates the Affleck-Dine leptogenesis in heavy 
gravitino/moduli scenario, which solves the moduli- induced gravitino problem while producing a 
correct amount of baryon asymmetry and relic dark matter density. The model involves two singlet 
flat directions stabilized by radiative corrections associated with supersymmetry breaking, one 
direction that generates the Higgs fi and B parameters, and the other direction that generates the 
scale of spontaneous lepton number violation. The dark matter is provided by the lightest flatino 
which might be identified as the axino if the model is assumed to have a U{1)pq symmetry to solve 
the strong CP problem. We derive the conditions for the model to satisfy various cosmological 
constraints coming from the Big-Bang nucleosynthesis and the dark matter abundance. 
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I. INTRODUCTION 



Low energy supersymmetry (SUSY) is one of the most plausible candidates for physics 
beyond the Standard Model (SM) at the TeV scale In the context of supergravity, 
SUSY can be spontaneously broken in a hidden sector while giving a vanishing cosmological 
constant. In this framework, the gravitino mass is related to the scale of hidden sector 
SUSY breaking as Msb ~ ^/m^pMpu where Mpi = 2.4 x 10^^ GeV is the reduced Planck 
mass. In low energy SUSY scenario, the hidden sector SUSY breaking is transmitted to 
the supersymmetric standard model (SSM) to induce soft terms providing sparticle masses 
of 0{1) TeV. Then, for certain range of ^3/2, gravitinos produced in the early Universe 
decay after the Big-Bang nucleosynthesis, which would destroy the successful prediction 
of the light element abundances [2|]. This cosmological difficulty might be avoided if the 
gravitino is relatively heavy, e.g. 7/13/2 > (9(10) TeV, so that the decay occurs before the 
nucleosynthesis. 

Recent progress in string flux compactification js^ has provided a SUSY breaking scheme 



in which the hidden sector with a SUSY breaking scale Msb ~ \frn^if2Mpi is naturally 

sequestered from the SSM sector. String fluxes generically produce a warped throat in 

compact internal space, and then the SUSY breaking sector is stabilized at the IR end of 

throat. On the other hand, the high scale gauge coupling unification suggests that the SSM 

sector is located at the UV end, and thus is sequestered from the SUSY breaking at the 

IR end of throat Ij, . In such case, the effective contact interactions between the SUSY 

il 

breaking hidden sector field and the SSM fields are so suppressed [6[ that the conventional 
gravity mediation [3] which would give a soft mass of 0(7/13/2) does not contribute to the 
SSM soft masses. As a result, the SSM soft masses can be much lighter than the gravitino 
mass, and 7713/2 = 0(10) TeV is a natural outcome if one takes the SSM soft parameters to 
be near the weak scale. Moreover, this framework can stabilize moduli either by fiuxes or 
by non-perturbative effects 0, and the resuhing moduli masses are typically much heavier 
than the gravitino mass, ~ 871^7113/2 for moduli stabilized by non-perturbative effects 
and ^ 87r^7n3/2 for moduli stabilized by fluxes [8|]. 

Nonetheless^the above scenario of heavy gravitino/moduli is not yet free from cosmolog- 
ical difficulty [9]. It has been noticed that there still arises a problem, the moduli-induced 
gravitino problem, due to that the gravitinos from moduli decays produce too many neutrali- 
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nos, which would overdose the Universe if the neutrahno is the stable lightest supersymmet- 
ric particle (LSP). As another difficulty, in heavy gravitino scenario it is not straightforward 
at all to get the Higgs ^ and B parameters having a weak scale size, which would require 
a severe fine tuning for the correct electroweak symmetry breaking. Recently, an attrac- 
tive solution to these problems has been proposed in involving a singlet fiat direction 
stabilized by radiative corrections associated with SUSY breaking. Once stabilized by ra- 
diative effects, the singlet fiaton S gets a loop-suppressed SUSY breaking F-component, 
F^/S ~ 7713/2/87?^, and then one can arrange the model to generate fi B /S, which 
would result in a weak scale size of /i and B when 171^/2 = 0(10) TeV. Such a scheme can 
solve the moduli-induced gravitino problem also. Now the LSP of the model is the fermionic 
component of S, the fiatino, which is much lighter than the MSSM neutralino. Then the 
neutralinos produced by moduli/gravitino decays are not stable anymore, but decay into 
light fiatinos with a relic mass density not overdosing the Universe. It is also natural to 
introduce a U{1)pq symmetry spontaneously broken by the vacuum value of S, and then the 
jseudoscalar component of S corresponds to the QCD axion solving the strong CP problem 



ll| , and the fiatino can be identified as the axino which has been proposed as a viable dark 



matter candidate in 
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On the other hand, late time thermal inflation [13| is a natural consequence of generic 
fiaton model, although its possibility was not explored in [l^. Thermal inflation also solves 
the moduli-induced gravitino problem by diluting away the coherent oscillation of moduli. 
A potential difficulty of thermal inflation is that any primordial baryon asymmetry is also 
washed away by the late time entropy production, so one needs a baryogenesis mechanism 
working after thermal inflation is over*. Recently, an interesting model of thermal inflation 
incorporating the Affieck-Dine (AD) leptogenesis [isj] has been proposed in the context of 
;he conventional gravity-mediated SUSY breaking scenario giving a weak scale size of ^3/2 



In gauge mediated SUSY breaking scenario, the baryon asymmetry produced by the Affleck-Dine mecha- 
nism can survive the dilution by thermal inflation if the initial amplitude of the AD flat direction is large 



enough [14[. In such case, the curvature of the potential at a field value larger than the messenger scale is 
determined by the small gravitino mass, thus the oscillation of the flat direction begins much later than 
the case of /2 ^ 'Tigoft , which would allow the generated baryon asymmetry large enough to survive the 
late time dilution. In our case of heavy gravitino scenario, this is not possible, so we need a baryogenesis 
after thermal inflation is over. 
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In this paper, we wish to examine the possibihty of thermal inflation incorporating a 
similar AD leptogenesis in heavy gravitino/moduli scenario realized in the framework of 
sequestered SUSY breaking. As we will see, a heavy gravitino mass requires additional 
conditions for successful thermal inflation and AD leptogenesis, e.g. the model should involve 
two singlet flat directions both of which are stabilized by radiative effects, one flat direction 
that generates the Higgs fi and B parameters, and the other flat direction that generates the 
scale of spontaneous lepton number violation. We present a viable model satisfying various 
cosmological constraints coming from the Big-Bang nucleosynthesis and the dark matter 
abundance. 

This paper is organized as follows. In section 2, we briefly discuss the moduli-induced 
gravitino problem and also how to generate the weak scale size of fi and B with a radiatively 
stabilized singlet flaton in heavy gravitino scenario. In section 3, we present a speciflc model 
of thermal inflation and AD leptogenesis in heavy gravitino scenario, and briefly analyze the 
phenomenological aspects of the model. A detailed discussion of the cosmological aspects of 
the model is provided in section 4, and section 5 is the conclusion. 



II. LOW ENERGY SUSY WITH HEAVY GRAVITINO 



A. Moduli-induced gravitino problem 

A gravitino mass much heavier than the SSM soft mass mgoft naturally appears in models 
with sequestered SUSY breaking. A SUSY breaking at the tip of throat in string flux com- 
pactification provides an attractive framework for sequestered SUSY breaking 0, ^. Such 
string compactiflcations include moduli that are stabilized by flux or non-perturbative effect. 
The resulting moduli masses are comparable to 87r^m3/2 for non-perturbative stabilization, 
and even much heavier than 87r^m3/2 for flux stabilization [§]. When SUSY breaking effects 
are taken into account, those moduli develop a nonzero F-component given by 



3/2 





where and denote the modulus vacuum value and the modulus mass, respectively. 
Since > 87r^m3/2, the resulting moduli-mediated soft masses do not dominate over the 



4 



anomaly-mediated soft masses of (9(7713/2/87?^)''' [8|, |22|, |23|. Thus, even when one takes 
into account the moduh-mediated soft masses of 0{F'^ /(p), it still holds true the order of 
magnitude relation mgoft ~ 7773/2/87?^, and the gravitino generically has a heavy mass of 
(9(87r^) Te"V when Trigoft is assumed to have a size of 0{1) TeV. 

Moduli are expected to have a misalignment of 0{Mpi) in the early Universe, and there- 
fore their coherent oscillations would soon dominate the energy density of the Universe. 
Since moduli masses are much heavier than (9(10) TeV in sequestered SUSY breaking sce- 
nario, their decays take place well before the nucleosynthesis. Nevertheless, moduli can 
cause a cosmological problem unless the branching ratio of the moduli decay to gravitinos 
is highly suppressed. Indeed, non-thermal gravitino production by heavy moduli decay can 
lead to a severe over- abundance of dark matter in the Universe 9|. 

For a modulus (f) with ^ '773/2, the number density of gravitinos produced by the 
modulus decay is given by 



"-3/2 _ 3 '^(ft^GG 



(2) 



S 2 r^^MSSM '^(/i 

where s is the entropy density, T^^^q and T^^mssm denote the decay width to the gravitino 
pair and to the MSSM particles, respectively, and is the reheating temperature for the 
modulus decays: 

(;^)'''v^.o.<-|^)''X^)-'^(^)--, (3, 

where g*{T) is the number of relativistic degrees of freedom at T, and = C^rn}^/ Mp^ is 
the total modulus decay width-!-. Gravitinos produced by modulus decay do not interact 
with others, and promptly decay into the LSP. The decay temperature of the gravitino is 
given by 

Tm - i^^S'^^W, ~ i ^Tt^, (4) 



.7r^5f*(T3/2)/ ^ ' 

with F3/2 ~ 10^^7773^2/^!'/ being the total decay width of the gravitino. Since 777,^ ^ 7773/2, 
T3/2 is much lower than T^. Under the assumption of R-parity conservation, if the LSP 



t In fact, moduli- mediation comparable to the anomaly mediation is phenomenologically desirable since the 
pure anomaly mediation [l9| suffers from the tachyonic slepton problem in the minimal supersymmetric 
SM ^ 



20 



21| 



^ For (/) "whose vacuum value determines the gauge coupling constants, it decays mainly into the gauge 
bosons and gauginos, and then — 0(10^^). 
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annihilation is not efficient, each gravitino will produce a stable LSP, yielding the LSP relic 
abundance given by 

S S ' S Vl-4X 10-5r<^^MSSM; ' ^ ^ 

where denotes the LSP mass, 

a = (^] (^) {^^\ {^^\ (6) 
VlO^V ^gATiv VlO^GevAlOSGeVy ' ^' 

and pcr/s — 1.9 x 10^^ GeV is the ratio of the critical density to the entropy density in the 
present Universe. 

Since Px/ Per in the present Universe cannot exceed about 0.25, the branching ratio of the 
modulus decay into the gravitino pair is bounded as 

r^^GG ^ , , ,,-5 flQ-'V" (9.iTt,W'' flO'GeV\ /lO^GeVV^' 



r^-^MSSM 



In view of that > 87r^m3/2 and 1713/2 = C(10) TeV in sequestered SUSY breaking scenario, 
r^^GG should be highly suppressed if x corresponds to the MSSM neutralino having a mass 
of O(IO^) GeV. However, it has been noticed that there is no suppression by 171^/2/712^ for 
the modulus decay into the helicity ±1/2 components of the gravitino [9]. As a consequence, 
even in the limit ^ 717^/2, the branching ratio is simply given by 



2 X 10-^ (8) 

n 

where A^^^ = 12 denotes the number of gauge bosons in the MSSM [9,]. This branching ratio 
exceeds the above bound by several orders of magnitudes when x is identified as the MSSM 
neutralino, and this is the moduli-induced gravitino problem. 

For a more careful analysis, one needs to include the effects of LSP annihilation. Including 
such effect, we find 

S S 4 \'K'^gAT3/2) J {(r^nnVTel)T3/2Mpi' 

where (aann'^^rci) is the thermal average of the annihilation cross section times the relative 
velocity of x- If the branching ratio of GG exceeds the bound ([7]), in order not to 
overdose the Universe, the annihilation is required to be as efficient as 

1.3x10-4/ 10 \ fWGeVy^^ 



On the other hand, even the Wino LSP has an annihilation cross section (cTann'^^rci) 
10~^ GeV~^, and thus it appears to be difficult to avoid the moduli-induced gravitino problem 
through the annihilation mechanism^. 



B. Higgs and Bfi term 

In heavy gravitino scenario, the MSSM Higgs B parameter generically receives a contri- 
bution of 0{m^/2) from SUGRA effect. Unless cancelled by other contribution, such a large 
B would make it difficult to realize the correct electroweak symmetry breaking. An attrac- 
tive way to avoid this difficulty is to generate fi and i? by a vacuum value of flat direction 



which is stabilized by radiative corrections associated with SUSY breaking 10|, l26|. Here 



we briefly discuss such a scheme using the model of [10| as an example. 

In addition to the canonical Kahler potential of the involved fields, the model includes 
the following Kahler mixing and the Yukawa interactions of the fiaton field S: 

AAnt = J d^9KS*S' + J £9^yH^H^Hd + ys^SS'j +h.c., (11) 

which are allowed by the global U{1) symmetry with the charge assignment: q{S) = q{S') = 
q{Hu,d) = 1, q{^) = —2. Then the direction along S ^ with S = S" = is fiat in 
the supersymmetric limit. After integrating out the heavy fields S and S' under a large 
background value of S, one obtains the effective theory of S described by 

A^eff = j d^e i^fS*S + (^K,f,^H^H, + h.c)j I , (12) 

with Yg^ = Ys{Q = \S\) and k^s = nyH/us- Here Q is the renormalization scale, and Ys is 
the running wave function of S. The potential is then determined by the running soft mass 

V = mUQ=\S\)\S\', (13) 



In fact, in other case that both and m^/o are of 0(10) TeV, the resulting Wino LSP abundance can be 



small enough not to overdose the Universe [2J,|25|. However, in this case, the modulus with ~ ™iZ2 



generically has F"^ /(f) TO3/2, and then the modulus-mediated sfermion masses can be of 0(10) TeV 25|. 
Here, we are focusing a different setup giving ^ 871^7713/2 (or heavier) with F'^/cj) ^ 7713/2/871^ (or 
smaller), which is motivated by sequestered SUSY breaking realized in string flux compactification. In 
our case, even the Wino LSPs produced by moduli/gravitino decays overdose the Universe by about two 
orders of magnitudes. 



7 



where m| = —F^ F'^*didj In Ys for denoting generic SUSY-breaking F-components in the 
model. The SUSY breaking F^ include first of all the F-component of the chiral compensator 
superfield 

c = Co + e^F^, 

as well as F^ and the moduli F-component F'^. For 5*0 = (S) > 771^/2, the equation of 
motion for F^ reads 

^ ~ -F'dj\nYs{Q = \S\). (14) 

Do 

It is quite possible that 1^5 is moduli-dependent at tree level, and then there will be a 
contribution of 0{F'^/(f)) to F^ /S. On the other hand, due to the super- Weyl invariance in 
the compensator formulation of SUGRA, the C-dependence of Ys arises only through the 
RG running. As a result, the contribution from F'^ is loop- suppressed, giving a contribution 
of 0{F^/8tt^Co) to F^Sq. In our case. 



F^ F^ "^3/2 ^ "^3/2 

Go rriff) ovr"^ 

and thus 

Now the fi and Bfi terms arise from the effective Higgs bilinear operator in the effective 
Kahler potential (|T2l) : 



where the loop-suppressed contributions from F'~^ are included in C(m3/2/87r^), together 
with the contributions from F'^'. Therefore, fi and B generated by a radiatively stabilized 
flat direction can be naturally of 0(m3/2/87r^), so be the weak scale when m3/2 = C(10) TeV. 



C. Axino dark matter and thermal inflation 

In heavy gravitino scenario, the addition of the singlet fiaton S has important implications 
for cosmology. Using the relation Yg^ = Ys{Q = \S\) and the stationary condition m'g{Q = 
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5*0) — for 5*0 > "^3/2 in the effective theory (fT2l) . one can find that the radial flaton and 
the flatino^ acquire SUSY breaking masses 

,2 



2 ^ dmliQ) 
- d\nQ 



= O 

-So 



"^3/2 



(Stt 



2^3 
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where 

and ^5 is the soft A-parameter associated with the Yukawa couphng ys in ( fTTI) . The angular 
flaton (= axion) as remains massless before the explicit breaking of U{l)s is taken into 
account. As having a light mass suppressed by a loop factor relative to the MSSM sparticle 
masses, the fiatino is the LSP and becomes a dark matter of the Universe under the usual 



assumption of R-parity conservation 12|, |27|]. Such a light dark matter is cosmologically 



favorable since the mo dull- induced gravitino problem would be considerably alleviated as 
can be seen in ([71). This point has been used in [l^ to solve the moduli- induced gravitino 
problem by assuming m^^ ~ 10^ MeV under the additional assumption that there is no 
thermal inflation triggered by the flaton S. 

However, since the flaton couples to the thermal bath in the early Universe through the 
Yukawa interaction which is responsible for its stabilization, it is a more plausible possibility 
that S is hold at the origin, 5* = 0, by its thermal mass in the early Universe. In that case. 



the Universe experiences a short period of thermal inflation (l3l . Il6l ] driven by the flaton 
potential energy Vq ~ itl'^^Sq around the origin. After thermal inflation, S would start to 
roll down towards its true minimum 5*0 and oscillate around the true minimum with an 
amplitude of 0{Sq). As it produces a tremendous amount of entropy, thermal inflation 
practically dilutes away all the unwanted relics, including the coherent oscillation of heavy 
moduli causing the moduli-induced gravitino problem. 

Thermal inflation is a natural consequence of the flaton model providing a weak scale size 
of fi and B in heavy gravitino scenario, and immediately solves the moduli-induced gravitino 
problem. However, there are certain cosmological constraints which require a considerable 



^ This fiatino can be called also the axino which would be the right name if the global symmetry U{l)s is 



good enough to be a U{1)pq symmetry solving the strong CP problem [ll| 
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extension of the model. After thermal inflation, the Universe is reheated by the decays of the 
oscillating radial flaton as- If the axion as is stable, the axion energy density produced by 
the decays of as is bounded by the Big-Bang nucleosynthesis, which requires that as decays 
dominantly into the SM particles 



28l | . In the model under consideration, to generate a weak 
scale size of -B, it is designed that S couples to the operator HuHd in the low energy effective 
action (fT2|) through the combination S* / S at tree level. With this structure of the model, 
the coupling of as to HuHd is cancelled at tree level, so its strength is loop-suppressed. As 
the decays of as to the SM particles are mediated by the coupling to HuHd, this results in 
that as decays dominantly to the axions, so an overproduction of axions which would be in 



conflict with the Big-Bang nucleosynthesis 28 1. 



Another difficulty to be overcome is that thermal inflation dilutes any pre-existing baryon 
asymmetry. Thus one needs to introduce a baryogenesis mechanism which works after 
thermal inflation is over. Recently, it has been noted that the observed baryon asymmetry 
can be generated via the Affieck-Dine (AD) leptogenesis after thermal inflation. The AD 
mechanism [isl can be implemented by a flat direction carrying nonzero lepton-number. 
A particularly interesting candidate for such a flat direction is the MSSM LHu as it can 
have a large nonzero value during the period of yU = 0, then rolls back to the origin after 
a nonzero /i is induced by the flaton vacuum value. Here, L denotes the lepton doublet 
superfield, and we do not specify the generation structure for simplicity. The dynamics of 
LHu is determined also by the dim = 5 neutrino mass operator: 

AW. = (18) 



which might be generated by the seesaw mechanism 29|]. Then, in order for the AD leptoge 



nesis to work, the A- type soft parameter associated with the neutrino mass operator should 



satisfy [la, ll8| 



\A,\^ < 6(m| + m^^ + |/i|2), (19) 

where rrii denote the soft scalar masses. Otherwise, LHu is trapped at a meta stable min- 
imum with LHu 7^ even after /i 7^ is generated. In heavy gravitino scenario, this is a 
nontrivial requirement as F'-^ of the SUGRA compensator generically gives a contribution 
of C(m3/2) to A.. 

As we will see in the next section, all these difficulties of thermal inflation in heavy 
gravitino scenario can be naturally solved by introducing an additional singlet flat direction 
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X. Thermal inflation can successfully incorporate the AD leptogenesis if X is stabilized 
also by radiative effects associated with SUSY breaking, and its vacuum value gives the 
heavy right-handed neutrino mass inducing the neutrino mass operator AWy via the seesaw 
mechanism. With this additional flaton X, the overproduction of axions is naturally avoided 
since now both fiatons, S and X, have an unsuppressed coupling to HuHj^, with which the 
fiatons decay dominantly into the SM particles. 



III. THE MODEL 



To resolve the difficulties noticed in the previous section, we introduce an additional 
singlet flaton X together with the right-handed neutrino N generating the neutrino mass 
operator AWi, via the seesaw mechanism, and an extra matter pair \t', \E''^ providing a Yukawa 
coupling to stabilize the original flaton S. The relevant part of the model is given by 



int 



d'^eKS*s' + I d^e 
+ I dH 



+ h.c. 



(20) 



where we have imposed two global U{1) symmetries, U{l)s and U{l)x, with the following 
charge assignments 

U{l)s : (2,2,0,-2,0,-1,2,-2), 
U{l)x : (0,0,2,-2,-1,0,2,0) 

for the fields S, S', X, S, N, L, HuHd and respectively. Obviously both X and 5* corre- 
spond to flat directions when S, S", N, \E' and \E''^ are all frozen at the origin. Here the extra 
matter fields \E', \1'^ can be either gauge-singlet or gauge- charged. In case of gauge-charged 
^>'^, our model can be considered as a simple generalization of some models discussed in 



30l |. incorporating thermal inflation and AD leptogenesis. 



Assuming that both X and S get large vacuum values, one can integrate out the heavy 
S, S", X, \E' and \E'^. The resulting effective theory can be written as 



eff — -^MSSM + / d'^6 
+ 



F|"X*X + + ( k^HuHd + h.c. 



X 



a -A,, h h.c. 

Z X 



(21) 
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where k = nyH/ux and Ay = y%/y'N- Here the effective wave function coefficient of ip = S, X 
IS given by Yf = Y^{Q = \ip\), where the running wavefunction coefficient Y^{Q) can be 
computed from the underlying theory (12U|) . 



A. Flaton stabilization 



For the underlying theory given by ( l20l) . the mixing between the two ffatons, X and S, in 
the effective potential is negligible. The ffaton potential is again induced by radiative effects 
associated with SUSY breaking, thus can be written in terms of the running soft masses: 

V ^ Yl K(Q = \^\M'^ (22) 

<fi=X,S 

where 

mJ(Q) = -F'F'*djdj\nY^{Q). (23) 

Here we consider the case that m'^{Q) is driven to be negative at certain low scales by the 
associated Yukawa interaction. Then the ffaton field is stabilized a.t (f = (fo satisfying 

mliQ = ^o) ^ 0, (24) 

and one can find the ffaton F-term is given by 

— ~ -F'di\nY^{Q = \ip\), (25) 
f 

where denote the SUSY breaking F-components in the model, including F^ ^ ^3/2 
and the moduli F-component ~ m^i^/m^ < 7713/2/871^. Again, for both 5* and X, the 
resulting F'^/y^o is of the order of the MSSM soft mass TTisoft ~ 1^3/2 /^t^"^ as desired. 

Around the minimum of potential, the radial ffaton a^p and the ffatino ip^p acquire SUSY 
breaking masses 



dmliQ) 



2 ""'"in 



ding 

2 

PK 

IGtt 



= o 



m 



2 

soft 



^7r2 



O (!^) , (26) 



where 7/^^ are the Yukawa couplings, A^ij are the soft A-parameters in the canonical basis. 
Here we are using the parametrization 
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for if = X,S, and also the fact that the mixing between X and S is neghgible. The 
axion component would remain massless unless the associated global U{1) symmetry is 
broken by higher dimensional operators or by non-perturbative effect. Because the flatinos 
have masses of (9(msoft/8vr^), the lightest flatino will constitute the dark matter under the 
assumption of R-parity conservation. 

B. Higgs fi and B, and the effective potential of LH^ 

In the model in consideration, the Higgs fi and B parameters are induced through the cou- 
pling of HuHd to the flaton field combination S* / X in the effective lagrangian fl2T]) . Includ- 
ing the potentially possible contribution from the moduli F-component ~ vn^^i^jm^ < 
"^3/2/877^, the resulting /x and B are given by 

^ = {(^) +^("^-*)}' 

B = — + OKoft), (27) 
^0 

where the contributions from F"^ are in 0{msoft)- One then finds both fi and B have a 
desirable size, i.e. 0{msoit) when kSo ~ Xq, which is rather a natural possibility. 

In the presence of the neutrino mass operator AW,^ = X^LHuLHu/2X , the effective 
potential for the LH^ flat direction is given by 

Vlh. = ^K + -k + H')l^l'+(^^ + --)+%^|T, (28) 

where i is the fleld variable parameterizing the LHu flat direction, and A^y is the soft A- 
parameter for AW^, which is given by 

Au = -T^ + 0{m,oit)- (29) 

During the period before a nonzero fi is induced by the flaton vacuum value, if m|+m|^^ < 0, 
which will be assumed in the following, i is stabilized at 

£0 ~ V|A,Xo/A,|. (30) 

If the condition f|T9l) is satisfied, which is easily done in our model, i is correctly rolling back 
to the origin after implementing the AD leptogenesis. Note that it is crucial for the AD 
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leptogenesis to work that Aj, is of the order of mgoft ~ ?7^'3/2/87r^, which is achieved in our 
model by generating the scale of lepton number violation by a radiatively stabilized vacuum 
value Xq. 



C. Flaton decay 

In this subsection, we examine the decay of the radial fiaton a^p {ip = S,T). Since 
nia^ — (!?(msoft/47r), these radial flatons are kinematically forbidden to decay into the MSSM 
sparticles, the gravitino, or the moduli. The axion components can acquire a small mass 
from non-perturbative effect or higher dimensional operator breaking U{1)^, but they will 
be taken to be almost massless in the following. In the low energy effective theory, the 
decays of into axions or fiatinos are mediated by the interactions 

^^^^ ^ ^ (^^Qy) (Q^Qy) + ^V-^ r!^^" cr^ip,pip<p + h.c, (31) 
2v2v?o 2V2(^o 

where A^^ is determined by the running A-parameters, A^jj, associated with the Yukawa 
couplings, 7/<^jj, responsible for the flaton stabilization: 

(32) 



Q='Po 



where i,j and k,l run over the fields having a nonzero Yukawa coupling with the flaton 
(p — X,S. One then finds that generically A^^ has a value of 0(10"^) or less. Then the 
decay rates are estimated as 

1 ml^ 
A?, mfi, m„ 

where we have taken kSq ~ Xq. 

The radial fiatons decay also to the SM particles, mainly through the effective interactions 

Ant = - I) + l/^H^^T - By^HlH^) + h.c., (34) 

which would induce a mass mixing between and the neutral Higgs fields after the elec- 
troweak symmetry breaking. The fiaton-Higgs mixing then allows a^p to decay directly into 
the SM fermions, which is dominated by the bottom quark channel: 



4n\ m\J \ mlj \mlJ 
14 



where b, h, and A are the bottom quark, the hghtest neutral CP even Higgs boson, and the 
neutral CP odd Higgs boson, respectively. 

After thermal inflation is over, the Universe is dominated by the energy density of coher- 
ently oscillating radial flaton a^. If decays dominantly to axions, it would be in conflict 
with the Big-Bang nucleosynthesis [28]. With the above results, the branching ratio to the 
decay into axions is given by 

^a^~^a^a^ 1 ( ^<t^\ ( ^'h\ (na\ 

T— ~ 48 \W> ' 

which can be easily smaller than O{10~^) to satisfy the bound from the Big-Bang nucle- 
osynthesis. 

Meanwhile, assuming Xq > So, the following interaction between two flatinos is induced 
by the X dependence of Yg^: 

Aatino = i^t^ i^xcr'^'ipsidf.as) +h.c., (37) 



where = {Xdx InYl''^). Then the heavier flatino ipi decays into the lighter flatino ipo plus 
an axion with the decay rate 

- ^y^-^J -x[- (38) 

Since A^ is induced at higher than two-loop level, its magnitude is suppressed as |A^| < 
(9(l/(87r^)^), so that ipi decays with a long lifetime while producing an axion together with 
the axino dark matter. Cosmological implication of such a late decay of the heavier flatino 
will be discussed in the next section. 



IV. COSMOLOGY 



The singlet flaton flelds in our model allow a natural realization of thermal inflation and 
baryogenesis, and provide light flatinos one of which is the LSP. We assume, as the boundary 
condition of our cosmology, that radiation is dominant before the commencement of moduli 
coherent oscillation, and the temperature was high enough to hold all flelds except moduli 
near the origin. As temperature drops down, the flaton fleld (f starts to roll down towards 
its true minimum ip = ipo a.t the critical temperature T^: 

T^ = ^ iv = X,S), (39) 
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' (p renormalized at 



3l| 



where = —rn?^{Q ~ mgoft) > is the soft scalar mass squared of 
Q ~ msoft, and [3^ determines how strongly ip couples to thermal bath 

Pi = lT.\y^^^\"- (40) 

ij 

Before the Higgs fi term is induced by the flaton vacuum value, £ parameterizing the LHu 
flat direction is also thermally trapped at the origin, and its critical temperature is given by 

T, = ^, (41) 

where mj = — (m| + m\j^)/2 > 0, and jSe is given by 

= l{j^\yuj\' + ^J2^^^(^^9l) + {L^H^), (42) 

ij a 

where ($) denotes the quadratic Casimir invariant of $. The pattern of thermal inflation 
and relevant cosmo logical contents depend on in what order the fields X, S and i roll down 
to the minimum. Here we consider the case that the underlying theory (1201) leads to 



Ts <Te < Tx, (43) 

which is a natural possibility. 

Meanwhile, the vanishing cosmological constant at the true vacuum {(f) = ipo, {£) = 
implies a nonzero potential energy near the origin: 

f=X,S <fi=X,S 

with 

^ (45) 

which has a value of O{10~^). Hence, during the epoch of the thermal confinement, Vq plays 
the role of vacuum energy. If necessary to be definite, we fix the vacuum values of the fiaton 
fields as 

Xo ~ 10^^ GeV, 5*0 ~ 10^^ GeV, (46) 

where we have taken into account that neutrino masses in the range of 10~^ eV are obtained 
for Xq ~ Aj^ X 10^^ GeV, and is of O(msoft) for kSo ~ Xq. 
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A. Thermal inflation 



When Hubble parameter H becomes comparable to the mass of the modulus at a 
time t — t^, the modulus starts its coherent oscillation with Planck scale amplitude, and 
the epoch of moduli domination begins. 

Thermal inflation begins aX t = t\ when the modulus energy density becomes compa- 
rable to Vo, and thus one finds 



This epoch of thermal inflation continues at least until X becomes unstable when the tem- 
perature drops to Tx Sit t = tx- 

The thermal history after tx crucially depends on the ratio between Vx and Vs- If 
Vx ^ Vs, the Universe is dominated by the non-rclativistic X particles. In this case, if 
Vs becomes dominant before S rolls down to its true minimum, a second epoch of thermal 
inflation is driven by Vs before or after i becomes unstable. Otherwise, matter domination 
by X particles would continue until the particles decay and reheat the Universe for the Big- 
Bang nucleosynthesis. On the other hand, if ^ Vs, thermal inflation is driven essentially 
by Vs, continuing until S becomes unstable, and the Universe is reheated by the decay of S. 

1. Vx:>Vs 

Since the commencement of the oscillation of X, its energy density px dominates the 
Universe as matter. As the Universe expands, px becomes comparable to V5 at i = ^2- 
The total decay width of the radial flaton (cp — X, S) can be written as 




(47) 



where pr is the energy density of radiation, and the temperature at ti is 




(48) 




(49) 



where has a value of (9(10 ^) because 



0{ml^j87r^). For £, one finds 




(50) 
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with Ci = 0(10-^) and 4 ~ ./\AJQK\ being the vacuum value of £ in the absence of the 
jj. term. For H ^V^, the contribution of px to radiation aXt — t2 leads to** 

TtM - i^oM)))' (^v.) ~ o.iy^r.M,, . (51) 

which implies 

For > 10^^ GeV, Tp^{t2) is lower than and thus the second thermal inflation can take 
place after I becomes unstable. Meanwhile, in case of Tp^(i2) < T^, the temperature due to 
the radiation contribution of at t — t2 is 

and therefore, by using V{1 — Qi) — V{1 — Iq) ~ m^^o A* = 0) we find 

TtAt^) ^ Ml\cJ\Cx ) \mxj ' 

-JT V^J iT^J iT^J feJ V^J • ^ ^ 

For < 10^^ GeV with Ps ~ ^(1)) the partial decay oi £ at t — t2 thus contributes 
dominantly to radiation with temperature larger than T5. This implies that, after £ decouples 
from thermal bath, a second epoch of thermal inflation begins at i = ^2 with a background 
temperature 

T2 = T,^{t2). (55) 

The thermal inflation ends as 5" rolls away from the origin, and then the decay of S eventually 
reheats the Universe for the Big-Bang nucleosynthesis^^ 

** We assume that the energy loss of the osciUating field due to the parametric resonance is small. 

The contribution to radiation from the late decay of X is negligible for our choice of the flaton vacuum 
values. 
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The number of e-folds for each epoch of thermal inflation is estimated as 



A^i ^ In 



T 



X 



No 



ilnf4^ 



~ 5.5 + 2 



6.8 + In 
1 



ax 



In 



10-1/ VlQi^GeV 



10=^GeV\^ flO^GeV\ 



rhx J 



J 



as 
10-1 



Tf 

2 



10-iy' ^10-3j \ ml )\XJ\ Xo ) 



(56) 



These e-folds, which are additive to the standard primordial inflation, do not interfere di- 
rectly with CMB or cosmic 21cm fluctuations"'""'", but may have an observable effect on the 
primordial density perturbation though it may be difficult to disentangle from uncertainties 
in the model of primordial inflation. 

The entropy productions from the decays of ip and i provide the dilution factors 



Ax ~ 
As ~ 



Tx 
Vs 



V: 



X 



PxiU 



Pe{ti 



Tx 



(57) 
(58) 
(59) 



where is the decay temperature of the flaton field S after thermal infiation. Therefore, 
the total dilution is given by 



A 



tot 



Ax A, Ac 



Vx V,Vl'^T,T\Mli 



T'xT, 



rpSrp4 

10-11 Cs 



1/2 



lOSo 

"x7 



m 



X 



moirLc 



103GeV\'/' , , 
) (60) 



ms 



1 /2 

where we have used ~ (F^Mp;) ' though we will redefine it later in a more definite way. 

Thermal infiation should dilute enough the abundance of moduli to avoid dark matter 
over-production described in Section [III The bound on the abundance at the time of the 



W For primordial inflation at V^/'^ ^ IQi" GeV and assuming radiation domination after inflation, if the 
reheating temperature of thermal inflation is Td = 0(1) GeV, the observable Universe leaves the horizon 
around 40 e-folds before the end of inflation. The observed CMB covers only around 6 e-folds after the 
horizon exit of the observable Universe 



around 9 additional e-folds 



33|. 



32l |. and the potentially observable cosmic 21cm fluctuations cover 
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decay of moduli is 



< 0.25-1^^ f ^) ^ 5.0 X 10---!^ ( '-^^) . (61) 



S ~ 'T^^^^rriy^ V S / present ^ ,j>-*GG \ 

For modulus particles produced before the thermal inflation, the late time abundance is 
estimated as 

^ ~ (^) ±- ~ 10-- (]^^\ (^^\ . (62) 
s J Atot \ J VAtot/ ' 

If produced at the end of the first thermal inflation, the modulus abundance is given by 
m^MliH{txY 1 1 



whereas, for those produced at the end of the second thermal inflation, we find 



m^Ml,H{tsY 1 VsT, 



d 



T| ml As mlMl, 



All of these contributions to the moduli abundance are far below the safe level of fl6ip . 



2. Vx<Vs 



In this case, we have only single epoch of thermal inflation driven essentially by Vs-, but 
it is extended by the radiation contribution of X and ^ when those fields decouple from 
thermal bath. The e-folding of the thermal inflation is given as 



iV ~ In 



Tx 



ax 



12 + - In 
3 



ax 



3 



Vc 







VxTxMpi VeTeMpi 



lO^GeV^ 



'4T/1/2 



T-14T/1/2 



,10-V V 10-4 J \ ml J 
The dilution factor from the decay of X is obtained as 



ms J 



lO^GeV 

m^ 



Ax 



Tx 



X 



pxite 



TeYVxTxMpi 



Tx) n v^/' 



1/2 



(65) 



(66) 
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while those from ^ and 5* have the same forms as (158|) and (l59l) . respectively. Hence the 
total dilution factor is 

. Vx V,T,TxMl, a\C, Cx .3 .4 ( rn,rn\ \ f So\ fMpA ^''^ 



.10- 

and thus, compared to the bound in (1^ . the abundance of moduli becomes totally negligible 
in this case too. 

After thermal inflation, the Universe is reheated by the decay of the flaton that ends the 
last thermal inflation. Since there is no unique decay temperature, we instead define the 
flaton decay temperature by 

Pr(Td) = ^rjM|„ (68) 

or 

^(7*(Td)Tl = PsM(Td) = ^r^^^sMr^Ml,, (69) 

where Fg-^^sM is the decay width to the SM particles, and ip is either X or S. This corre- 
sponds to a time 

^ (70) 

and at the moment we have 

p^(Td) ^ ^rj, (71) 

and the entropy increases after by a factor Sf/Sd — 2 where and Si are the entropy at 
t = td and a late time, respectively [isl . 



B. Baryogenesis 

For Tx > Ti > Ts, the flaton field X decouples first from thermal bath and settles to 
its true minimum X = Xq. The scalar potential implementing the AD leptogenesis is then 
determined by the effective theory fl^T]) . which involves 

Cad = I d'eYfSS*+(^J d'9k^H^H, + j d'O^X^^^^^^ + h.c)j , (72) 
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where Yg^ = Ys{Q = \S\). To analyze the dynamics induced by m'j^{Q) + mfj^i^Q) < at 
Q ~ msoft, we parameterize the associated flat directions'^ as 



Hd = (0, hd 



(73) 



with the D-term constraint for the SU{2)l gauge symmetry 

D2 = \hu\'-\hd\'-\l\' = 0. 
The relevant part of the potential is then written as 

Vad = Vs + ^glDl + ml\S\^ + ml\l\^ + mjjJKl 

(\ g* 
+ -A.XJ^hl - Bfi—Khd + 
\Z Ao 



(74) 



m 



2 \h |2 



c.c. 



S* 

\yfhu + fi—K 

^0 



to 


S* 


+ 


^0 



(75) 



with fi given by 



S 



+ 0(msoft; 



(76) 



where the loop-suppressed contribution from F'^ is included in (9(msoft)- Note that B, Ay, 
and fi have values of (9(msoft) and are nearly independent oi ip = X,S for \lp\ ^ mgoft since 
the equation of motion for F"^ is given by (!25l) . 

Initially, all fields are held at the origin by their finite temperature potential. As the 
temperature drops down, one of the unstable directions, S oi C. = Ih^, will roll away from 
the origin. We assume Ihu rolls away first, i.e. Ti < T5. Then the potential term ^A^X^Phl^ 
fixes the phase of Ihu, while |A^//;.^p and jAjy/^/iup stabilize its magnitude. The Ihu field 
may partially reheat the thermal bath and so prolong the thermal inflation, but eventually 
S will also roll away from the origin, ending thermal inflation. As S rolls away, the term 
B(iS*huhd/ will force hd to become non-zero. This provides temporarily a large mass to 
quark and lepton directions and constrains those directions to zero, shielding the dynamics 
from the dangerous non-MSSM vacuum in the direction associated quark and lepton 181.134 1. 



Then, BfiS*huhd/Xo fixes the phase of S*huhd. As S nears its minimum, the cross term 



We set quark and lepton directions to be zero since they are expected to be held at the origin throughout 
the dynamics. This will become clear from the subsequent argument. 
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from \Xi,Phu + piS*hd/ Xq\'^ rotates the phase of Ihu generating a lepton asymmetry, and at 
the same time \ jlS*hu/ Xq\'^ gives an extra contribution to the mass squares of Ihu and /i^^d, 
bringing them back in towards the origin. Thus, we have a type of the Affleck-Dine (AD) 
leptogenesis. Preheating then damps the amphtude of the Ihu and h^h^ fields keeping them 
in the lepton preserving region near the origin 17|]. The Ihu and huhd fields then decay, at 
a temperature in the MSSM sector above the electroweak scale, and their lepton number is 
converted to baryon number by sphaleron processes. Finally, the fiatons 5* and X decay, 
diluting the baryon density to the value required by observations, ub/s 10^^°. 

One may think that, if X decays later than S, the baryon asymmetry generated due to 
the dynamics of S may be significantly diluted. However, the dilution is possible only when 
Vx ^ Vs with Xq different from 5*0 by more than several order, which is very unnatural in 
our model. 

Comparing to the model considered in 18|], the Higgs /i term in our model has a linear 
dependence on the triggering field S rather than quadratic: 

S* 

= At^, (77) 

where fi does not depend on 5 for l^l ^ J^j-soft- This may weaken the strength of torque 
responsible for the angular momentum and preheating, but would only make a difference 
of a factor of 0{1). We therefore expect that the AD baryogenesis would work well in our 
model too. 



C. Dark matter 

The model (I2T]) provides the lightest flatino as the dark matter under the assumption of 
R-parity conservation since flatinos ip^ [if = X,S) have small masses of (9(msoft/87r^). In 
addition, the axions can also become a dark matter of the Universe if they are light enough 
to be stable. These dark matter components can be cold, warm or even hot, depending on 
their masses and how they are produced. In this subsection, we will derive cosmological 
constraints on the dark matter. 
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1. Cold dark matter 



As a cold dark matter, the lightest flatino is dominantly produced by the decay of the 
corresponding flaton if it eventually reheats the Universe for the Big-Bang nucleosynthesis. 
Another source is the decay of heavier sparticles in thermal bath, which is dominated by 
the decay of the lightest MSSM sparticle. In addition, axion misalignment and strings also 
contribute to the energy density of cold dark matter. 

a. Flatinos produced by the flaton decay: For the flaton (f that ends the thermal infla- 
tion, its decay reheats the Universe. The late time flatino abundance from the flaton decay 
is determined by [l8| 



2r, 



dt' a^{t')p^{t') 



2.2Td r, 



where the decay temperature is roughly given by 



Td ~ IGeV 



100 



lo- 



rn,, 



103GeV J 



y / 10^^GeV \ 



(78) 



1 1/4 



(79) 



which follows from fl69|) . Using fl33|) and fl78l) , one can find the current abundance of flatino 
dark matter ipo 



a 



5.6 X 10* 



8 / "^^0 A ^ 



IGeVJ s 



3.6 



100 



1/2 



10 



-1 



IGeV) VlGeV/ V / 



2 /lGeV\ /10iiGeV\ 



(80) 



where we have used that one flatino dark matter is produced per each heavier flatino ipi 



if ipy, = ipi, and that the total decay width is given by r,^ 
requirement Q^,^ < ^^cdm — 0.2 translates into 



+SM- 



Therefore, the 



< l.SGeV 



100 



IGeV J \ X 



10-1 



lOi^GeV 



(81) 



which is at the lower end of its expected range (12611 . 

b. Flatinos produced by the decay of thermally generated MSSM sparticles: The flaton 
coupling to HuHd in the effective Kahler potential fl2T|) induces 



Ant = Xn^Hu,dHd,u^s + i^Hu4Hd,u(^''df,'ilJs - -^Hu^Hd^ui'x + h.c, (82) 
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where = (Sdslnk) = C(l/87r^), A'^ = kSq/Xq = 0{1), and Hu,d are the Higgsino fields. 
The above interactions lead the lightest MSSM sparticle x to decay into fiatinos with the 
decay rate |l2[ | 

C.j, TTvt , ^ 

where C^g = (9(10^^C^^) for /i of the soft mass scale and Sq ~ lO^-'^Xo. Here C^^ ~ 1 may 
contain a factor of m|/m^, and we have neglected the masses of decay products. 
The thermal bath generates x with the number density 



oo 



= ^ I dk == , (84) 

and they subsequently decay into fiatinos, dominantly into ipx- The late time fiatino abun- 
dance is thus estimated as 



whose numerical solution can be found in llSl 



s 

where F^{x) is approximated as 



, 5.3x^ for X <^ 1 
F^x) ~ { , (87) 

5.4 for X > 1 



with 



and ~ 2m^/21 being the temperature at which the fiatino production rate is maximized. 
Using (!83l) and (!86|) . the current abundance of fiatino dark matter is obtained 

^ " ^ -^'^-(^w^) (^) (^) (^)^^*'^'''-' 

where we have used T^, ^ Fo-^^^sm, and that the number density of ipo is the same as that of 
ipx even when ipx = '^i- Therefore, for x <^ 1, one finds 
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which does not exceed ^cdm if the flaton decay temperature satisfies 

103 )\m^)\ m,A„ / V 1013C 



(91) 



6.7 V5'*(^d) 

Since Td is expected to be a few GeV, the above cosmological bound is well satisfied for 
= O{10^) GeV. Note that the fiaton decay temperature can be less than or similar to 
the freeze-out temperature of which is about [32]. Therefore, in order to avoid 

direct production of x from the fiaton decay, we may require mo-^ < 2m^, which is quite 
plausible in our model. 

Flatinos will also be produced by the decay of x after they freeze out. However, the 
standard Big-Bang neutralino freeze-out abundance is good match to the dark matter abun- 
dance, our freeze-out abundance of x will typically be less than the standard abundance, 
and m^^ -C m^, therefore the fiatino abundance generated after the freeze-out should be 
safe. 

c. Axion misalignment: The cold dark matter contributions of misalignment and ax- 
ionic strings from QCD-axion are well known. Thus we consider only the cases of non-QCD 
axion under the assumption that the mass of the axion appears due to a tree-level symme- 
try breaking term which is small enough not to disturb the radiative stabilization of the 
associated fiaton field. In this case, the mass of axion is not connected with its coupling 
constant. 

The energy density of an axion misalignment when oscillation commences is 

pa - 2 AT^' > 

where a = a^, is the axion mass, is the misalignment angle, and M is the vacuum 
degeneracy. The current energy density is 



Pa f ao \ g*{Tosc) ( g*s{To) V''^ / Pa 



Pi V'^osc/ 5'*(To) \g*s(Tosc) J \Pt 

"'g'J\T,sc) f gM) y^'e' mlv' 



g* (To) V 9*s (Tosc) ; Af^ {3HoscMpif^ To ' 
Therefore, Qa < ^cdm requires 



(93) 
(94) 



< 6.0 X lO-^eV — — - — , 95 

\g*{Tosc)J V J 

where we have used 3ifosc = and g^{Tosc) = g*s{Tosc)- 
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Thermal inflation can dilute those contributions if the reheating temperature is lower 
than the temperature at which the mass of axion becomes comparable to the expansion 
rate. We do not analyze the dilution in this paper, instead we refer the reader to |3| for 
the case. 



d. Axionic strings: The axion produced by the strings is 42| 



A- 



ToscMpi 

where A = [lor In (/i/iJosc)]- The current energy density is 



(96) 



Pa 4 /vr^ ,^ 
7. = 3^U'*^ °' 



1/4 



9*{To) J 3HolcM%^To' 
Therefore, in order for Qa not to exceed f^cDM, it is required 



(97) 



rua < 2.5 X 10-^eV 



50AV 100 flO^^GeVy 



(98) 



where we have again used SHqsc = ^a- 

The mass bound fl98|) implies that tree level symmetry breaking should be very small 
For instance, the higher dimensional superpotential term 



AW oc 



n-3 ' 
PI 



[n > 3) 



(99) 



provides the axion a mass 



ml ~ \A^o\ 



M 



PI 



n— 3 



(100) 



where A ~ ^3/2 (or mgoft). In this case, the requirement ( l98i) is satisfied for 
\n{ml/\A^o\) 



n >3 + 



PI 



> 7.4 In 
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(101) 



2. Warm/hot dark matter 



The LSP ipQ and axion can be warm or hot when they are produced by the decay of 
the next to LSP (NLSP) tpi (for LSP and axion) and the radial flaton a^^ (for axion). Th e 



constraint on hot dark matter comes from CMBR and structure formation 



35 



36 



33,1381 
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Currently allowed hot dark matter fractional contribution to the present critical density is 



nHDM < 10- 



39|. The constraint may be more stringent as suggested from the analysis 



of the early re- ionization of the Universe at high redshift 40|. Taking into account the 



recent analysis of WMAP 5-year data 



4l| . the allowed warm/hot dark matter fractional 



contribution to the present critical density is likely to be 

^^WHDM < 10-'. (102) 

We will take this as the upper bound on the fractional energy density of our warm/hot dark 
matter. 

a. Hot axion from the flaton decay: Axions produced by the flaton decay have a current 
momentum 

Pa = ^, (103) 

where a is the scale factor at the time they were created and Oq is the scale factor now. 
Whereas, the current momentum of an axion produced at id is 



«o 2 Si/'gli\Td)Td 2 - > v^ys^^^^^^ 



so it may be relativistic now. The current number density spectrum is given by 



dpa V«o/ H Hpl 

which may provide an observational test of our model in the future. The energy density of 
the axions is 

_ 9.{Td)g%\T) r.^ 



(106) 



Therefore, assuming that the hot axions are still relativistic now, their current energy density 
is estimated as^^ 



100 \^ r,. 



^ 2 X 10^^ ( (107) 



The energy density of thermally produced axions is fit, ^ ma/lSleV, hence it will be subdominant 42l ]. 
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b. Hot axion from the decay of the NLSP: Because of the coupling (137|) . axions can be 
produced also by the decay of the NLSP, i.e. the heavier flatino ipi, which is originated from 
the decays of flaton and x- Neglecting axion mass, the current momentum of the axion from 
the NLSP*** is given by 

p. = - ,108) 

ao 2m^i 

In particular, the current momentum of the axion produced at t'^ = T^^^^^ is 

Pd 



ao 2m^, 9l'siT'^)Tk ^^^^ 




^ 1.87 X 10-^eV -rjP^ 1-^1 (109) 

where is the background temperature when the NLSP decays 

T', - (^^.ra)"''(r,,^,„M„)^/^ - lOeV [-^^Y- (110) 

In this case, the axion would thus be ultra-relativistic now for any plausible axion mass. 

The current number density spectrum is 

p dnl"^ _ f aV 2p^, r^,^^o _ Wa ^ ^i^i^^oVl ^^^^^^ 
" dpa \aoJ H m^^ Hp'J^ 

which again may provide an observational test of our model in the future. The late time 

energy densities of the axions are 

g^{Td)gts{T) \ m^, P^ 



^hot 
ra 



PSM 



flaton \g.{T)g,s{TA)g%\K) J T-s-sm 



„hot 
ra 



(112) 
(113) 



where {n^i/s)t'^ is given by ( 186|) . and their current energy densities are 

»l°'L.,„ - 4 X 10-= f-^V (114) 



The energy densities in ( I107p . (11141) and (11151) are typically well below the bound of (I102p 
in our model. 



We ignore the effect of NLSP's momentum, since NLSPs from the decay of Batons and x '^iU be highly 
non-relativistic during most of its life. 
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c. Warm/hot LSP from the decay of NLSP: Eqs. f llOSp - (11 111) are applicable in this 
case too. Note that, although it is heavy and non-relativistic now, the LSP may be warm or 
even hot unless the NLSP (ipi) and LSP (ipo) are highly degenerate, and may have observable 
astrophysical effects. The late time energy density of the LSP is simply given by 

^Vo ~ — —^ipi, (116) 



where Q,^^ is the would-be energy density of the NLSP ipi if it did not decay. Therefore, 
from dHU), ([89]), ([90]), ffT02D and ffTTSD . Q^^ < ^whom requires for NLSPs from fiaton 




where F,^ ~ r„ ^sm has been used. 



V. CONCLUSION 



Heavy gravitino with 7713/2 = 0(10) TeV is a generic prediction of sequestered SUSY 
breaking scenario. String flux compactification provides a natural setup for sequestered 
SUSY breaking, i.e. SUSY breaking at the IR end of warped throat, and also can stabilize 
all moduli by fluxes or nonperturbative effects. The resulting moduli masses are much 
heavier than the gravitino mass, ~ 871^7713/2 or even heavier, which might be considered 
as an attractive feature in view of the cosmological moduli problem. However such heavy 
moduli still cause a cosmological difficulty, the moduli-induced gravitino problem, producing 
too many LSPs which would overdose the Universe if the LSP is given by the MSSM 
neutralino. Another potential difficulty of sequestered SUSY breaking scenario is that it is 
not straightforward to get a weak scale size of the Higgs fi and B parameters. 

An attractive way to get a weak scale size of /i and B in heavy gravitino scenario is to 
generate them by a singlet flat direction which is stabilized by radiative effects associated 
with SUSY breaking. Unless one assumes an unusual type of initial condition, such a flat 
direction generically triggers a late thermal inflation. Thermal inflation would immediately 

30 



solve the moduli-induced gravitino problem, however it requires a baryogenesis mechanism 
to work after thermal inflation is over. 

In this paper, we have presented a model of thermal inflation in heavy gravitino sce- 
nario, which successfully incorporates the Affleck-Dine leptogenesis while producing a cor- 
rect amount of relic dark matter density. The model involves two singlet flat directions 
stabilized by radiative effects, one direction that triggers thermal inflation and generates a 
weak scale size of fi and B, and the other direction that generates the scale of spontaneous 
lepton number violation. The dark matter is provided by the lightest flatino which might 
be identified as the axino if the model is assumed to have a U{1)pq symmetry to solve the 
strong CP problem. The collider signal of the model highly depends on the flaton vacuum 
values (ip) {if = X, S) which determine the rate of the decay of the lightest sparticle in the 
MSSM to the lighter flatinos. If (if) are well above the intermediate scale ~ 10^° GeV, the 
collider phenomenology of the model would be almost the same as the conventional MSSM. 
For {(f) ~ 10^'^ GeV or lower, the model can give a distinct signal associated with the decay 
of the lightest sparticle in the MSSM and this will be the subject of future work 43]. 
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